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CONSTRUCTION  OF  LARGE  REACTION  MECHANISMS 

INTRODUCTION 

Researchers  are  currently  studying  combustion  reactions  in  an  effort  to  gain  new  insights  into  fuel 
combustion  efficiency  and  fuel  safety.  The  understanding  of  combustion  systems  is  often  greatly  aided 
through  the  use  of  computational  models.  However,  the  complexity  of  multispecies  systems  makes  the 
modeling  task  formidable.  The  methodology  presented  here  demonstrates  the  logical  formulation  of  a 
reaction  mechanism  for  a  multispecies  combustion  system. 

For  the  combustion  modeler  and  the  chemical  modeler  in  general,  the  formulation  of  a  complex 
mechanism  consists  of  two  major  steps: 

1  Specification  of  the  molecules  and  molecular  fragments  to  be  included  in  the  system. 

2.  Selection  of  the  chemical  reactions  to  be  considered  and  included  in  the  proposed  mechanism. 
Subsequent  steps  to  complete  the  overall  modeling  process  include: 

3.  Assignment  of  the  rate  constants  to  the  selected  reactions. 

4.  Solution  of  the  time  dependent  coupled  differential  equations  of  the  model's  full  mechanism. 

5.  Comparison  of  the  output  from  the  model  with  experimental  data. 

6  Performance  of  sensitivity  analysis  on  the  parametric  input  to  the  model. 

The  six  steps  are  usually  performed  in  the  order  listed,  but  not  necessarily.  This  report  is  con¬ 
cerned  mostly  with  Step  2,  the  selection  of  the  chemical  reactions  in  the  mechanism.  We  will  assume 
that  the  species  involved  have  been  previously  specified  (Step  1)  and  we  will  use  a  sample  species  set 
for  illustrative  purposes. 

Frequently  chemical  kinetic  mechanisms  are  the  product  of  a  researcher’s  expertise  and  feelings. 
However,  it  can  be  very  difficult  for  the  modeler  to  consider  each  of  the  possible  reactions  when  con¬ 
structing  a  reaction  set  from  a  given  set  of  reactive  species.  Even  a  small  species  set  can  give  rise  to  an 
extremely  large  reaction  scheme.  Typically  combustion  systems  can  have  more  than  30  species.  If  the 
modeler  does  not  or  cannot  consider  all  possible  reactions  involving  the  given  species,  then  it  is  possi¬ 
ble  that  potentially  important  reactions  will  not  be  considered. 

Usually  there  is  some  agreement  among  modelers  as  to  the  rate  constants  to  assign  to  various 
reactions  (Step  3).  Many  models  are  attacked  because  of  rate  constants  that  are  possibly  incorrect  by 
one  or  two  orders  of  magnitude  (typical  units  of  cm’  molecule  1  s  ').  However,  one  quickly  realizes 
that  omitting  a  potentially  important  reaction  underestimates  that  rate  constant  by  13  or  14  orders  of 
magnitude.  (Omission  sets  a  rate  constant  of  zero  versus  a  typical  bimolecular  rate  constant  value  of 
I01'  to  I01'1  cm1  molecule  1  s  '.)  It  is  for  this  reason  that  it  is  quite  important  to  consider  (not  neces¬ 
sarily  include,  but  consider)  all  possible  reactions.  Even  the  most  expert  kineticist  can  easily  err  when 
it  comes  to  mechanisms  of  50  or  more  species.  It  should  be  mentioned  parenthetically  here  that  the 
difficulty  in  the  eventual  solution  of  the  coupled  differential  equations  (Step  4)  is  in  the  number  of 
species,  not  the  number  of  reactions,  because  the  number  of  coupled  differential  equations  is  dictated 
by  the  number  of  species.  A  model  of  900  reactions  of  55  species  is  tractable. 
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If  a  systematic  procedure,  or  algorithm,  which  allows  for  the  consideration  of  all  possible  reactions 
can  be  followed  when  constructing  a  reaction  set,  some  of  the  uncertainty  associated  with  the  proposed 
mechanism  can  be  eliminated.  The  organization  of  the  mechanism  is  also  enhanced,  and  the  time 
required  to  construct  the  complete  reaction  set  is  decreased.  Therefore,  in  an  effort  to  reduce  the  arbi¬ 
trariness  inherent  in  Step  2,  we  present  in  this  report  a  convenient  algorithm  by  which  modelers  can 
quickly  and  automatically  construct  a  complete  reaction  set.  We  assume  for  this  analysis  that  all  the 
relevant  species  which  play  a  significant  role  in  the  mechanism  have  been  specified. 

SAMPLE  SYSTEM  AND  NOTATION 

In  order  to  understand  the  use  of  the  proposed  algorithm,  it  will  be  discussed  here  in  the  context 
of  an  example  system.  Specifically,  we  formulated  and  used  our  algorithm  as  part  of  our  effort  model 
understanding  the  low-temperature  n-butane  oxidation  phenomenon.  Some  of  the  criteria  explained  in 
the  following  section  are  specific  !o  the  n-butane  oxidation  system;  however,  most  considerations  are 
generally  valid. 

The  set  of  reactive  species  used  as  a  sample  in  this  report  is  listed  in  Table  1.  The  column  titled 
"Species"  gives  each  species  an  eight  character  name  (four  computer  words  on  our  16-bit  machine). 
The  name  is  for  the  convenience  of  the  users;  no  exclusion  criterion  is  based  on  the  name  of  the 
species.  Columns  2  through  4,  NC,  NH,  and  NO  designate  the  number  of  carbons,  hydrogens,  and 
oxygens  in  the  original  molecule.  Remember  that  the  system  of  interest  here  is  a  C,  H,  O  system.  The 
algorithm  can  be  easily  modified  to  include  other  atomic  species.  Column  5  entitled  NZ  gives  a  number 
to  each  species.  This  number  is  mostly  for  internal  accounting  in  the  program.  It  is  shown  in  the  next 
section  how  selective  numbering  can  aid  in  generating  an  additional  criterion.  Columns  6  and  7  are  the 
AHf(298)  in  kcal/mole  and  S(298)  cal/(deg  mole),  respectively,  of  the  input  molecules.  These 
numbers  can  be  found  in  the  literature  [1]  or  estimated  using  conventional  means  [2]. 

It  will  be  helpful  to  understand  the  following  discussion  of  the  algorithm  and  the  program  itself, 
to  introduce  some  notation.  The  set  of  stable  molecules  and  molecular  fragments  (radicals)  shall  be 
called  R  There  are  N  species  in  R  including  a  special  species  M  which  is  explained  later  in  this  report. 
Let  RU).  R(J),  R(K),  and  R(L)  represent  the  I,  J,  K,  and  L  species  within  the  set  of  the  N  chosen 
species.  Indices  1,  J,  K  and  L  range  independently  over  the  entire  set  of  N  species.  There  is  nothing  to 
prohibit  these  indices  from  being  equai.  Throughout  the  text,  and  the  program,  I  and  J  will  index  reac¬ 
tants  and  K  and  L  will  index  products.  Thus,  ail  possibilities  of  two  reactants  forming  two  products  are: 

R  ( I )  +  R  ( J)  —  R(K)  +  R(L).  (1) 

We  can  represent  Eq,  (1)  using  matrix  notation.  The  reactants  portion  of  Eq.  (1)  are  represented 
by  the  N  x  N  matrix  whose  entries  consist  of  all  possible  combinations  of  R (I)  and  R(J).  Since  there 
are  N  possible  R(I)  and  independently  N  possible  R ( J) ,  there  are  N2  possible  species  combinations  of 
R(l)  +  R(J) 

As  a  simple  example,  consider  the  species  set  consisting  of  only  A,  B  and  C.  We  can  represent  all 
possible  combinations  of  reactants  involving  these  species  by  the  entries  in  the  following  symmetric 
matrix: 

AA  AB  AC 

BA  BB  BC  . 

CA  CB  CC 

Similarly,  this  matrix  also  represents  ail  possible  combinations  of  products  of  bimolecular  reactions 
involving  these  species.  The  special  species  M  is  included  in  the  reaction  set  so  that  addition  and  uni- 
molecular  decomposition  reactions  are  also  available. 
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Thus,  all  possible  reactions  in  the  form  of  Eq.  (1)  can  be  represented  as  a  mapping  from  one 
entry  of  the  reactant  matrix  onto  any  one  entry  of  the  product  matrix.  For  the  species  set  A,  B,  C,  the 
following  are  equivalent  for  the  set  R: 


AA 

AB 

AC 

AA 

AB 

AC 

BA 

BB 

BC 

— ►  BA 

BB 

BC 

CA 

CB 

CC 

CA 

CB 

CC 

There  are  N2  entries  in  each  matrix  and  thus  N2  x  N2  —  N4  total  reaction  possibilities.  Even  for 
a  relatively  small  set  of  species,  the  total  number  of  possible  reactions  is  formidable.  There  is  an  obvi¬ 
ous  need,  therefore,  to  formulate  criteria  to  eliminate  many  of  the  possible  reactions  quickly  and 
automatically. 


CRITERIA  IN  MECHANISM  GENERATION 

Appendix  A  lists  program  RXNGN,  a  FORTRAN  program  for  generation  of  a  reaction  mechan¬ 
ism.  The  many  criteria  appearing  in  the  program  RXNGN  will  be  discussed  in  this  section  but  not 
necessarily  in  the  order  that  they  appear  in  the  program.  Generally,  the  criteria  in  the  program  appear 
in  order  of  decreasing  exclusion.  The  sooner  a  reaction  fails  a  criterion  and  can  be  removed,  the  less 
criteria  must  be  applied  to  that  reaction.  The  development  of  criteria  in  the  text  will  be  ordered  more 
along  the  lines  of  a  chemical  deductive  process. 

One  of  the  most  important  criteria  for  the  reactions  to  satisfy  is  mass  balance.  The  input  informa¬ 
tion  on  the  number  of  carbons,  hydrogens,  and  oxygens  (columns  2,  3,  and  4  of  Table  1 )  is  used  for 
the  mass  balance  computation  (lines  83  to  89  of  the  program).  As  a  programming  note,  we  parentheti¬ 
cally  point  out  that  carbon  balance  is  checked  fust,  then  hydrogen  balance,  and  then  oxygen  balance.  If 
the  carbon  balance  fails,  there  is  no  need  to  check  the  hydrogen  or  oxygen  balance  (see  program). 

The  number  of  possible  reactions  that  may  be  eliminated  by  the  mass  balance  criterion  is  entirely 
dependent  on  the  specific  species  list,  and,  therefore,  cannot  be  denoted  generally  in  terms  of  the 
number  of  original  species,  N.  The  following  discussion  indicates  the  number  of  reactions  that  are 
eliminated  by  the  imposition  of  each  criterion,  but  the  number  indicated  does  not  exclude  those  already 
eliminated  by  the  mass  balance  criterion.  Note  that  the  number  of  reactions  that  may  be  eliminated  on 
the  basis  of  any  one  criterion  is  highly  dependent  on  the  order  in  which  the  criteria  are  imposed. 

The  next  criterion  to  consider  is  that  of  simple  permutation.  Since  it  does  not  matter,  for  exam¬ 
ple,  if  A  reacts  with  B,  or  B  with  A,  all  of  the  following  denote  the  same  reaction  chemically: 

A  +  B  —  C  +  D 

A  +  B  —  D  +  C 

B  +  A  —  C  +  D 

B  +  A  —  D  +  C. 

This  repetition  needs  to  be  eliminated  systematically.  To  illustrate  the  systematic  elimination,  we  use 
the  simple  matrix.  Label  the  upper  and  lower  triangular  region  of  each  symmetric  matrix  X|,  x2,  x3,  x4 
such  that: 


AA 

AB 

AC 

AA 

AB 

AC 

BA 

BB 

BC  — ► 

BA 

BB 

BC 

CA 

CB 

CC 

CA 

CB 

CC 

becomes 


h 


I 


I 
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AA 

,4.4 

*4 

BB 

— » 

BB 

X) 

*3 

CC 

CC 

Since  the  product  and  reactant  matrices  are  symmetric  about  the  main  diagonal,  X!  and  x2  of  the  reac¬ 
tant  matrix  and  x3  and  x4  of  the  product  matrix  contain  equivalent  entries.  To  eliminate  the  repetition, 
we  elimi”ate  all  the  entries  of  either  x(  or  x2  and  x3  or  x4.  If  we  arbitrarily  choose  to  eliminate  x2  and 
x4,  we  no*  have: 


AA 

AA 

BB 

— ► 

BB 

X| 

*3 

CC 

CC 

By  eliminating  the  entries  of  x2  and  x4,  we  eliminate  (N  -  1)  +  (N  -  2)  +  (N  -  3)  +  ...  +  1  — 
N(N  -  1) 

- - -  entries  from  each  matrix.  The  number  of  reaction  possibilities  we  have  eliminated  by 

imposing  this  criterion  is: 


The  number  of  entries  eliminated 
from  the  reactant  matrix. 

The  number  of  entries  remaining 
in  the  reactant  matrix. 


Total  number  of  entries! 
in  the  product  matrix.  J 

The  number  of  entries  eliminated 
from  the  product  matrix. 


which  in  terms  of  N  is: 

N(N  -  1) 


N2  +  N2  - 


N(N  -  1) 


n(n  -  i; 


N2(N  -  1) 


(3N  +  1) 


reaction  possibilities  eliminated.  The  number  of  entries  remaining  in  each  matrix  is  N2  - 
N(N  -  1 )/ 2  —  N(N  +  l)/2.  Since  there  are  (N)(N  +  1 )/ 2  entries  remaining  in  each  matrix,  and  since 
with  the  remaining  entries  all  mappings  are  possible,  there  are 

N(N  +  1)  N(N  +  1)  f  N(N  +  1) 12 


2  2 
possible  reactions  remaining  to  be  considered. 

Simple  permutations  are  eliminated  by  the  code  on  line  91  in  the  program  in  Appendix  A. 

Nonreactions  or  "identity"  reactions  (i.e.,  A  +  B—  A  +  BorA  +  A—  A  +  A,  etc.)  must  be 
eliminated.  Thus,  our  criterion  is  that  if  I  =  K  and  J  =  L,  or  I  =  L  and  J  *■  L,  or  1  —  K  and  J  —  K, 
or  I  =  L  and  J  —  K,  then  the  reaction  is  eliminated  from  consideration.  These  reaction  types  are  exact 
mappings  from  an  entry  in  the  reactant  matrix  onto  the  same  entry  in  the  product  matrix.  Since  there 
are  the  same  number  of  these  identity  reactions  as  there  are  remaining  entries  in  each  matrix,  by  elim¬ 
inating  these  identity  reactions  from  consideration,  we  eliminate  N(N  +  1 )/ 2  possible  reactions.  The 
number  of  remaining  possible  reactions  is  then 

N(N+1)  | 2  N(N  +  1)  <N  -  1 )  (N)  (N  +  1)(N  +  2) 
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The  next  applicable  criterion  is  no  reactions  with  certain  species  or,  if  on  the  basis  of  experimental 
evidence  it  can  be  determined  that  certain  species  appear  in  negligible  relative  concentrations  or  are 
unreactive,  then  one  may  want  to  consider  those  species  as  products  only,  or  eliminate  them  from  con¬ 
sideration.  If  a  certain  species  is  to  be  considered  as  a  product  only,  it  is  necessary  to  eliminate  all  pos¬ 
sible  reactions  involving  that  species  as  a  reactant  going  to  all  possible  products.  For  any  given  species 
there  are  exactly  N  remaining  entries  in  the  matrix  which  contain  that  species.  For  example,  if  we  are 
interested  in  species  A ,  we  find  that  for  the  case  where  N  =  3,  there  are  exactly  3  entries  which  contain 
at  least  one  A. 

AA 

BA  BB 

CA  CB  CC 


The  number  of  possible  product  entries  to  which  a  given  reactant  entry  can  go  is  the  number  of 
entries  on  the  product  side  minus  one.  We  subtract  one  from  the  number  of  entries  on  the  product 
side  to  avoid  double  counting  those  reactions  which  we  eliminated  previously  by  imposing  the  identity 
reaction  criterion.  For  example,  if  we  wish  to  eliminate  all  reactions  involving  species  A  as  a  reactant, 
then  we  do  not  want  to  count  the  following  reactions  as  being  eliminated  by  this  criterion: 

A  +  A  —  A  +  A 
B  +  A  —  B  +  A 
C  +  A  —  C  +  A. 

Thus,  removing  a  given  species  from  consideration  as  a  reactant  eliminates: 

Number  of  entries 
N  ■  —  1 

on  product  side 

_  N  .  NiN_Lii_  , 


N2  +  N  —  2 


possible  reactions. 


If  for  some  reason,  after  the  species  set  has  been  chosen,  we  wish  to  to  eliminate  all  possible  reac¬ 
tions  involving  a  given  species  as  both  products  and  reactants,  then  we  eliminate  the  possible  reactions 
involving: 


A  reactant  pair  having  at  going  products  minus 
least  one  of  the  pairs  as  ►  exact  mappings 
species  to  be  eliminated.  10 

Product  pairs  which  have 
The  remaining  going  at  least  one  of  the  pairs 
reactant  pairs.  (  being  the  species  to  be 
eliminated 


Or  in  terms  of  N  we  eliminate: 

N  - - - I  +  - - - 1N=»  (N'  —  1)  ■  N  possible  reactions. 

(For  an  example  of  this  criterion  see  lines  1 12  and  1 13  in  Appendix  A.) 

Because  radicals  usually  are  more  reactive  than  stable  species,  reactions  involving  two  radicals,  or 
one  radical  and  one  stable  species  as  reactants,  are  much  more  likely  to  occur  at  low  temperatures  than 
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reactions  involving  two  stable  species.  A  fifth  possible  criterion  is  to  eliminate  all  possible  reactions 
involving  two  stable  species  as  reactants.  This  criterion  is  applicable  to  our  low  temperature  oxidation 
of  n-butane,  but  is  not  applicable  in  many  other  cases.  Letting  NSP  represent  the  number  of  stable 
species  in  the  species  set,  we  find  that  the  number  of  stable  reactant  pairs  in  the  reactants  matrix  is  NSP 
(NSP+D/2.  We  now  wish  to  eliminate  all  reactant  pairs  consisting  of  two  stable  species  which  go  to 
all  possible  products  minus  exact  mappings.  In  our  simple  example,  let  us  assume  A  is  the  only  stable 
species. 


AA 

AA 

BA  BB 

BA  BB 

CA  CB  CC 

CA  CB  CC 

We  wish  to  eliminate  the  following  possible  reactions: 

A  +  A  —  B  +  A 
A  +  A  —  B  +  B 
A  +  A  —  C  +  A 
A  +  A  —  C  +  B 
A  +  A  — *  C  +  C. 

(Note:  A  +  A  —  A  +  A  was  eliminated  under  the  exact  mapping  criterion.)  Here  (with  N  =  3  and 
NSP  —  1)  we  eliminate  five  possible  reactions.  In  general  with  NSP  stable  species  and  N  total  species, 
we  eliminate 

NSP(NSP  +  1)  N(N  +  1)  _  . 

2  2 

"possible  reactions." 

This  criterion  is  applied  in  the  program  in  Appendix  A  by  line  104.  The  species  in  the  input  file 
were  carefully  grouped  (column  NZ)  so  that  all  stable  species  were  numbered  after  all  the  radicals,  and 
the  special  species  M.  In  Table  1  there  are  29  stable  species. 

It  is  useful  to  limit  the  reactions  also  on  an  enthalpy  and  free  energy  basis.  With  the  input 
enthalpy  (AW/T298),  column  6  of  Table  1)  and  entropy  S(298),  column  7  of  Table  1)  it  is  easy  to  cal¬ 
culate  A//#  and  AG*  (lines  131  to  133  of  program  RXNGN  in  Appendix  A): 

AHr-  AHP  -  AH, 

ASr  “  ASP  -  AS, 

AGr  “  AHr  -  TASr, 

where  T  represents  temperature,  and  the  subscripts  R,  p,  r  represent  reaction,  products,  and  reactants. 

One  can  then  set  limits  for  acceptable  AGR  and  AHr.  Ideally,  of  course,  AGr  should  be  less  than 
or  equal  to  zero.  However,  since  measurements  of  AHf(298)  and  S(298)  are  not  perfect,  and  since 
some  estimated  values  for  those  quantities  were  used,  it  is  advisable  to  set  the  cutoff  for  AGR  above 
the  theoretical  limiting  value  of  0.0. 

In  large  complex  chemical  systems,  such  as  combustion,  endothermic  reactions  occur.  In  our 
example  program  we  used  a  cutoff  value  for  AH*  of  35  kcal/mole  (line  134  of  program  RXNGN  in 
Appendix  A). 

The  next  two  criteria  we  used  involve  a  knowledge  of  the  structure  of  each  molecule.  In  a  car¬ 
bon,  hydrogen,  oxygen  system  there  are  five  major  types  of  bonds:  1)  carbon-hydrogen  (the  number  of 
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carbon-hydrogen  bonds  is  designated  NCH);  2)  carbon-carbon  (NCC);  3)  carbon-oxygen  (NC'O);  4) 
oxygen-oxygen  (NOO);  and  5)  oxygen-hydrogen  (NOH).  Oxygen-oxygen  bonds  are  important  in  low- 
temperature  9-butane  oxidation  because  of  the  prevalence  of  peroxides;  this  bond  type  might  be  unim¬ 
portant  in  other  systems.  Logically,  three  atomic  species  should  give  rise  to  six  bond  types.  The  sixth 
type  here  is  hydrogen-hydrogen,  and  since  that  occurs  in  only  one  molecule  in  our  system  hydrogen 
<H2),  NHH  is  not  included.  A  double  bond  would  count  as  two  bonds.  For  example,  NCC  of  ethene 
is  2.  Table  2  lists  the  numbers  of  each  of  these  five-bond  types  for  our  input  data  set.  The  first 
column  (NZ)  corresponds  to  the  NZ  column  of  Table  1.  (Note  that  Table  2  concatenated  to  Table  1  is 
essentially  the  input  file  to  program  RXNGN.) 

Using  the  number  of  each  type  of  bond  as  input  data,  there  are  two  apparent  restrictions  to 
impose.  First,  limit  the  total  bond  change  between  products  and  reactants;  and  second,  put  limits  on 
the  change  in  each  type  of  bond  between  products  and  reactants. 

It  must  be  kept  in  mind  that  the  purpose  is  to  generate  a  set  of  elementary  (one-step)  reactions. 
If  the  products  have,  for  example,  four  carbon-carbon  bonds  and  the  reactants  have  nine,  it  is  difficult 
to  see  how  one  elementary  chemical  step  can  break  five  carbon-carbon  bonds.  The  number  of  each 
type  of  bond  broken  is  subtracted  from  the  number  of  each  type  of  bond  formed  (lines  138  to  142  of 
program  RXNGN  in  Appendix  A).  Specific  limits  are  then  imposed  for  the  change  in  total  number  of 
bonds  (lines  148  and  152  to  154)  and  the  change  in  number  of  each  type  of  bond  (lines  157  to  161) 
(Appendix  A). 

A  special  species  M.  without  any  mass,  ( i.e. ,  the  number  of  carbons,  hydrogens,  and  oxygens  in  it 
are  all  zero)  is  used  to  allow  for  unimolecular  decompositions,  addition  reactions,  and  single  counting 
of  isomerization  reactions. 

Unimolecular  decomposition:  A  +  M  — *  B  +  C. 

Addition  Reaction:  A  +  B— •  C  +  M. 

Isomerization:  (if  A  and  B  are  isomers  then  all  of  the  following  satisfy  mass  balance): 

A  +  M  —  B  +  M 
A  +  B  —  B  +  B 
A  +  C  —  B  +  C  ,  etc. 

We  eliminate  all  isomerizations  except  for  the  first  reaction  listed.  So  that  the  program  is  more  general, 
M  does  not  have  to  be  a  specific  species  number  (NZ),  thus  when  numbering  the  species,  no  special 
number  is  reserved  for  M.  The  program  finds  what  the  species  number  of  M  is  (lines  117  to  122). 
The  species  number  of  M  is  not  used  for  the  unimolecular  decompositions  or  additions.  The  redun¬ 
dancy  of  isomerizations  is  eliminated  by  the  code  on  lines  124  and  125  of  Appendix  A. 

CRITERIA  APPLIED  TO  n-BUTANE  OXIDATION 

We  now  want  to  examine  the  results  of  applying  our  general  criteria  developed  in  the  first  section 
of  this  report  to  the  species  set  which  we  considered  in  modeling  the  low-temperature  oxidation  of  n- 
butane.  Our  species  set  consisted  of  56  species  (27  radicals,  28  stable  species,  and  the  special  species 
M,  see  Table  I).  As  we  noted  earlier  for  the  reaction  R (I)  +  R ( J )  — •  R(K)  +  R(L)  with  the  only  res¬ 
triction  being  that  1,  J,  K,  and  L  be  members  of  our  species  list  consisting  of  N  species,  there  originally 
existed  N4  or  in  this  case  (56)4  =  9,8348,496  possible  reactions  to  be  considered. 

With  56  species  both  the  reactants  matrix  and  the  products  matrix  were  originally  56  x  56  square 
symmetric  matrices  composed  of  3136  entries  each.  By  imposing  the  permutation  criterion,  (which 
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essentially  eliminates  all  entries  above  or  below  the  main  diagonal  from  each  matrix)  we  eliminated 
N(N  -  1 )/ 2  =  56(55)/2  =  1540  entries  from  both  the  products  and  reactants  matrix.  The  number  of 
reaction  possibilities  which  we  eliminated  was: 

N2— [3N  +  11  =  (56)2  •  —  |3(56)  +  11  =7,287,280. 

4  4 

The  number  of  reaction  possibilities  which  remained  was  (N(N  +  l))2  =  2,547,216. 

By  imposing  the  identity  reaction  criterion,  we  eliminated  N(N  +  l)/2  possible  reactions  or,  in 
this  case,  56(57)/2  =  1596  reactions.  Thus  by  imposing  the  first  two  criteria,  we  eliminated  7,290,416 
possible  reactions. 

For  our  particular  system  we  selected  three  stable  species  to  act  as  products  only.  As  shown  ear¬ 
lier,  removing  a  given  species  from  consideration  as  a  reactant  results  in  the  elimination  of 
N(N2  +  N  -  2 )/ 2  possible  reactions  or  in  this  case: 

(three  species)  ( 56) 562  +  56  —  2)/2  =  267,960 
possible  reactions  were  eliminated 


We  imposed  the  criterion  to  eliminate  all  possible  reactions  involving  two  stable  species  as  reac¬ 
tants.  Since  in  our  particular  species  set  there  were  28  stable  species,  by  applying  this  criterion  we  elim¬ 


inated: 


NSPINSP  +  I) 
2 


N(N  +  1) 


28  m) 
2 


56(57) 


-  1  =  647,975  possible  reactions. 


Before  taking  the  mass  balance  criterion  into  consideration  and  as  a  result  of  imposing  these  four  cri¬ 
teria.  our  reaction  set  consisted  of  1,628,145  possible  reactions.  In  our  example  the  mass  balance  cri¬ 
terion  was,  next  to  the  permutation  criterion,  the  most  restrictive.  Every  species  in  our  system  was 
composed  of  some  number  and  combination  of  II  atoms.  C  atoms,  and  O  atoms.  The  mass  balance  cri¬ 
terion  was  then,  if  one  of  the  following  was  not  true,  the  reaction  was  eliminated: 


C(K)  +  C(L)  -  C(J)  -  C(l)  =  0 
IKK)  -I-  H(L)  -  H(J)  -  H ( I )  =  0 
C) ( K )  +  OIL)  -  O(J)  -  0(1)  =  0. 


Using  only  three  criteria,  mass  balance  and  elimination  of  permutation  and  redundancy  reduced 
the  set  from  ( 56 )4  =  9,834,496  to  a  much  more  manageable  13,092. 

Imposition  of  the  thermodynamic  criteria  (limiting  value  for  AHR  and  AGR)  reduced  the  reaction 
set  by  3282  reactions.  The  bond  breaking  and  forming  criteria  reduced  the  set  by  2687  reactions.  By 
allowing  isomerizations  to  occur  only  with  the  special  species  M,  the  mechanism  shortened  by  another 
1080  reactions.  The  specific  case  criterion  of  no  bimolecular  reactions  of  stable  species  eliminated  539 
reactions.  The  additional  specific  checks  (see  listing  of  program  in  Appendix  A)  eliminated  2255  reac¬ 
tions  and  brought  the  set  down  3249  reactions.  Remember,  the  amount  each  criterion  reduces  the  total 
is  dependent  on  the  order  in  which  the  criteria  are  applied. 

Thirty-two  hundred  reactions  are  still  quite  a  lot,  but  few  enough  so  that  visual  inspection  is  a 
realistic  task.  We  examined  the  mechanism  with  the  aid  of  the  search  programs  which  are  described  in 
the  next  section.  Many  of  the  remaining  reactions  were  not  one-step  reactions.  In  fact,  many  reactions 
appeared  to  be  two-step  reactions  and  both  of  the  individual  steps  were  already  in  the  mechanism. 
Below  are  three  such  examples  where  part  (a)  is  an  obvious  two-step  mechanism  and  parts  (b)  and  (c) 
are  the  two  individual  reaction  steps  which  were  already  part  of  the  mechanism. 
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1)  a.  CHjCO  +  02  -  CHj02  +  CO 

b.  CHjC'O  +  M  —  CH3  +  CO 

c.  02  +  CHj  —  CH302  +  M 

2)  a.  CHjCO  +  n— C3H7  — »  C4III()  +  CO 
b  CHjCO  +  M  — ’  CHj  +  CO 

c.  CH3  +  n-C,H7  —  C4H,()  +  M 

3)  a.  CHjOO  +  i-C,H,  —  i-C3H702  +  CH3 
b  CH,00  +  M  —  CH3  +  02 

c.  0:  +  i-C3H7  —  i-C3H702  +  M 


The  output  of  program  RXNGN  is  not  reproduced  here  because  it  is  excessively  long;  however, 
the  outputs  of  the  two  search  programs  described  in  the  next  section  are  formatted  the  same.  (See 
Tables  3  and  4.)  For  each  reaction,  the  reactants  and  products  with  their  respective  species  numbers; 
the  AIIK;  the  AGK.  the  number  of  reactions  tried,  the  number  of  reactions  found;  the  net  change  of 
C  H.  CC,  CO,  OO,  and  OH  bonds,  and  the  total  net  bond  change,  are  tabulated.  The  species  numbers 
are  for  convenience  and  give  the  subsequent  search  programs  something  to  key  on.  The  number  tried 
is  the  number  of  reactions  that  pass  the  mass  balance,  redundancy  and  permutation  criteria.  The 
number  found  represents  those  reactions  that  have  passed  all  the  imposed  criteria.  The  bond  change 
numbers  represent  the  sum  of  that  type  of  bond  in  the  products  minus  the  sum  of  that  type  of  bond  in 
the  reactants.  (Note  that  these  numbers  can  be  negative.)  The  column  entitled  NET  is  the  sum  of  the 
absolute  values  of  each  of  the  preceding  five  columns. 

From  the  original  (56)4  =  9,834.496  possible  reactions,  we  were  able  to  construcl  a  reasonable 
reaction  set  by  formulating  a  logically  ordered  set  of  general  and  specific  rational  criteria,  and  imposing 
these  criteria  upon  all  of  the  possible  reactions.  We  are  thus  able,  to  a  great  extent,  to  remove  the  pro¬ 
cess  of  reaction  set  formulation  out  of  the  intuitive  domain,  imposing  the  formulated  rational  criteria 
allowed  us  to  systematically  arrive  at  a  reaction  set  consisting  of  about  850  reactions,  after  originally 
considering  over  nine  million  possible  reactions 

SEARCH  CAPABILITIES 

Three  programs  were  developed  to  assist  in  the  examination  of  the  mechanism.  These  three  pro¬ 
grams  look  at  the  file  created  by  the  program  RXNGN  (Appendix  A). 

The  capability  to  search  through  the  data  file  is  very  powerful.  For  example,  program  SELET 
(reproduced  in  Appendix  B)  searches  the  reaction  file  for  all  reactions  in  which  a  particular  input 
species  is  found  as  a  reactant,  or  as  a  product;  or  as  either  a  reactant  or  a  product.  In  this  way,  one  can 
quickly  see  the  formation  and  destruction  pathways  a  particular  species  takes  through  the  overall 
mechanism.  If  the  user  is,  for  example,  interested  specifically  in  methanol,  the  reaction  set  may  be 
searched  for  the  reaction  in  which  methanol  appears.  Alternatively,  one  could  search  for  those  reac¬ 
tions  that  create  hydroxy  radicals.  It  is,  of  course,  possible  that  from  examination  of  this  edited  output 
from  the  search  program,  the  chemist  might  feel  that  a  reaction  that  is  felt  to  be  important  intuitively  is 
missing.  This  means  that  either  (I)  a  key  species  has  been  omitted  from  (he  original  data  file  input  to 
RXNGN;  or  (2)  one  of  the  criteria  used  in  RXNGN  is  too  stringent  and  excludes  important  reactions. 
Chemical  knowledge  and  intuition  are  still  important  in  the  use  of  these  algorithms. 
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Table  3  shows  some  sample  output  of  program  SF.LET.  All  the  reactions  in  which  methanol 
(species  number  28)  occurred  as  a  reactant,  were  found. 

Program  SLRXT  (reproduced  in  Appendix  C)  permits  the  user  to  search  for  reactions  of  classes 
of  reactants.  For  example,  one  could  examine  all  the  reactions  that  occur  with  alkyl  radicals  by  input¬ 
ting  the  species  numbers  of  all  the  alkyl  radicals.  The  program  also  allows  complementary  reactants, 
that  is,  one  could  examine  all  the  reactions  that  occur  between  alkyl  radicals  and  alcohols  by  specifying 
all  the  alkyl  radicals  and  (hen  specifying  all  the  species  numbers  of  all  the  alcohols  as  complementary 
reactants.  Table  4  shows  this  situation  as  a  sample  output  of  SLXRT 

The  last  utility  program  is  entitled  RUST  (reproduced  in  Appendix  I))  This  simple  program 
allows  the  listing  of  the  reaction  file  in  a  shortened  version  so  that  there  is  room  (in  the  right  hand 
margin)  to  annotate  the  file.  The  other  two  programs  (SFLKT  and  SLRXT)  reproduce  the  input  file 
exactly  as  it  is,  merely  selecting  out  certain  reactions.  Program  RUST  is  merely  a  convenience  and  is 
included  as  Appendix  I)  for  completeness. 

CONCLUSION 

We  have  formulated  an  algorithm  to  systematically  produce  tractable  reaction  mechanisms  for 
computational  modeling  of  complex  chemical  kinetic  systems.  The  additional  searching  algorithms  pro¬ 
vide  a  very  powerful  observational  and  educational  tool  as  they  permit  the  user  to  study  thoroughly  the 
entire  reaction  mechanism 
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LIST  OF  SYMBOLS 
Description 


A,  B,  C,  D 

Designation  for  arbitrary  species 

1,  J 

Indices  of  reactant  species 

K,  L 

Indices  of  product  species 

M 

Special  species 

N 

Number  of  species  in  set  R 

NC 

Number  of  carbon  atoms 

NCC 

Number  of  carbon-carbon  bonds 

NCH 

Number  of  carbon-hydrogen  bonds 

NCO 

Number  of  carbon-oxygen  bonds 

NH 

Number  of  hydrogen  atoms 

NHH 

Number  of  hydrogen-hydrogen  bonds 

NO 

Number  of  oxygen  atoms 

NOO 

Number  of  oxygen-oxygen  bonds 

NOH 

Number  of  oxygen-hydrogen  bonds 

NSP 

Number  of  stable  species 

NZ 

Species  number  of  a  specific  species 

R 

The  set  of  all  input  species 

R(index) 

The  indexed  species 

S  ( 298 ) 

Entropy  at  298K 

T 

Temperature 

x,x2x,x4 

Areas  of  particular  matrix 

under  discussion 

sgr 

Free  energy  change  of  reaction 

AH,  (298) 

Heat  of  formation  at  298K 

Enthalpy  change  of  reaction 

Subscript  p 

Products 

Subscript  r 

Reactants 

Subscript  R 

Reaction 

Units 


cal/(degmole) 
degrees  K 


kcal/mole 

kcal/mole 

kcal/mole 
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Table  1  —  Species  and  Parameters 


SPECIES 

HC 

HH 

HO 

HZ 

DELTA  H 

S 

0* 

4 

4 

1 

24 

39  3 

38  3 

02 

4 

4 

2 

48 

0 . 0 

49  0 

H* 

4 

1 

4 

1 

32  .  1 

27  .  4 

OH* 

4 

1 

1 

2 

9  .  4 

43  .9 

HO  2* 

4 

1 

2 

3 

3 . 0 

54 . 4 

H20 

4 

2 

1 

49 

-37.8 

43  .  1 

CH4 

1 

4 

4 

30 

-17  9 

44  .  3 

CO 

1 

4 

1 

44 

-26.4 

47  .  3 

C02 

1 

4 

2 

40 

-94 . 0 

31.2 

C  H  30  * 

1 

3 

1 

S 

3.5 

55 . 2 

CH30H 

1 

4 

1 

28 

-48.  1 

37 . 3 

C  H  302* 

1 

3 

2 

6 

6  .  7 

65  3 

C2H4 

2 

4 

0 

29 

12.5 

52  .  4 

C9H5* 

a 

5 

4 

7 

at .  s 

SO  .  0 

C  2  HE 

2 

6 

4 

30 

-20 . 2 

54 . 9 

CH3C0* 

2 

3 

1 

8 

-5.4 

64  .  5 

CH3CH0 

2 

4 

1 

31 

-39.7 

63 . 2 

CH3CH20* 

2 

5 

1 

9 

-4.0 

65  .  3 

C2H500* 

2 

S 

2 

10 

-3.0 

76  4 

C3H6 

3 

6 

4 

32 

4  .  9 

63  .  8 

C  3  H  7* - M 

3 

7 

4 

11 

210 

68  3 

C3H7*-  I 

3 

7 

0 

12 

17  6 

66  7 

C2H3CH0 

3 

4 

1 

33 

-17.7 

69  1 

CH3C0CH3 

3 

6 

1 

34 

-51.7 

70  5 

C3H70-I* 

3 

7 

1 

13 

-12.6 

71  8 

C4H8TBTE 

4 

0 

4 

33 

-2  .  7 

70  9 

C4HB1BTE 

4 

8 

4 

36 

0 . 0 

73.6 

C4H9-H* 

4 

9 

4 

14 

16.8 

76.6 

C4H9-S* 

4 

9 

4 

13 

13  7 

75 . 2 

C4H10HBT 

4 

14 

4 

37 

-30  .  2 

74  .  1 

C4H80THF 

4 

a 

1 

38 

-44  .  0 

73  6 

C  4  M  BOH 

4 

14 

1 

39 

-63  .  6 

86  8 

C4H90* 

4 

9 

1 

16 

-14  0 

76  5 

C4H902*S 

4 

9 

2 

27 

-17.8 

93  5 

H 

4 

0 

4 

4 

0  .  0 

0  4 

H  2  02 

4 

2 

2 

4  1 

-32  .  6 

36  4 

HC  MO 

1 

2 

1 

42 

-26 . 0 

32  .  3 
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1 

1 
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17 

9  0 

53 . 7 

CH300H 
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4 

2 

43 

-31.3 

67  3 
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26 

34  .  3 
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-17.5 
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88 . 2 
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C3H700H 
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8 
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84  .  6 
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8 

2 

33 
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82  .  7 

C3H8 

3 

8 

4 

34 

-24  8 

64  .  3 

C4H823E8 

4 

8 

1 

36 

-23 . 9 

84  6 

H2 

4 

2 

4 

33 

4  0 

31  .  2 
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Appendix  A 
PROGRAM  RXNGN 


i 


PAGE  000  1  FTN  9:3a  fin  THU  ,  16  JULY,  1981 


0001  F  TN4  ,  L 

00)2  PROGRAM  RXNGN<3.99),  REV  800709  OOREH  IH0R1T2 

0003  C  VERSION  800709  FOR  H  6  H  L  E  T  7  PACKARD  MINICOMPUTER 

0004  C  IDEA  BUILT  UPON  MALT  SHhUB'S  REACTION  GENERATOR 

0005  C 

0006  C  PROGRAM  GENERATES  ALL  POSSIBLE  REACTIONS,  GIVEN  INPUT  SPECIES 

0007  C  (UP  TO  100),  SUBJECT  TO  VARIOUS  CONSTRAINTS  LISTED 

0008  C  THROUGHOUT  THE  COMMENTS  IN  THE  PROGRAM 

<10  0  9  C 

0010  C  CURRENTLY  DIMENSIONED  FOR  UP  TO  100  INPUT  SPECIES 

0011  DIMENSION  H ANE!  4  ,  100 ) , EN T L ( 1 0 0  ) , HZ!  10 0 > , EH TR( 1 0 0  ) 

0012  DIMENSION  M C H<  1  0 0  ) , N CC ( 1 0 0  )  ,  NCQ(  10 0  ) , HOO ( 1 00  ) ,  H 0H(  1 0 0 > 

0013  INTEGER  C<  100  )  , H<  1 0 0  ) , Q(  1  00  ) ,  I  P< 5  ) 

0014  OATA  I  I  ,  J J  ,TEMP70 , 0  ,  299  07 

00  13  C 

0016  C  FORMATS 
0017  C 

0018  100  FORMAT!  ’  ENTER  LU  OF  DATA  CASSETTE)  ') 

0019  101  FORMAT!  IX, 619) 

0020  102  FORMAT!  '  IER  ■  '  ,  18  ) 

0021  103  FORMAT!  12  ) 

00  22  104  FORMAT! 1 X ,  4  A  2 , 4  X ,  I1,3X,12,4X,  I1,4X,I2,2X,F103,2X,F10.3) 

0023  105  FORMAT! 2X, 4A2, 4N,  I1,3X,12,4X,  II, 4X,  I2.2X.F10  1.2X.F10  11 

0024  106  FORMAT!  '1  REACTION  GENERATOR  H-BUTAHE  COOL  FLAME' 7  77  ) 

0023  107  FORMAT! 2X, 4A2,  '(', 12, ')’. 3H  *  ,4A2,'(',I2,')'#4H  -)  , 

0026  +  4A2  ,  ' < '  ,  12  ,  '  )  ' , 3H  4  ,  4 A 2 . ' < '  ,  I  2 , '  > '  . 2X , F 10  1 . 2X , 

0027  +  F10  .  1  ,  IX  ,  15,  16 , 615  ) 

00  28  108  FORMAT!  '1  '  ,9X,  'REACTANTS'  ,22X»  'PRODUCTS'  ,18X,  '  0  E  L  H  '  ,  8X ,  'DELG', 

00  29  t  3X,  1 U  TR  D ’ , 2X,  ’ DFHD ' , 2X  , ’NCHD’ , IX, *  HOC  0  * . IX, ' NCOD ' . IX, 

0030  8  'NOOO', IX, 'N0HD',2X, 'NET'  ) 

0031  109  FORMAT!  ’  LEN  ■  110  ) 

0032  110  FORMAT! 7777'  TOTAL  REACTIONS  EXAMINED  «  ',110) 

0033  111  FORMAT! 613  ) 

0034  112  FORMAT!  1H1 ,7X, 'HZ'»6X,  'NCR'  ,6X,'HCC',6X,  'HCO'  . 6  X , 'HOO' ,6X,'N0H'  > 

0033  113  FORMAT!'  ENTER  LU  OF  OUTPUT  DEVICE') 

0036  114  FORMAT!'  RXNGN  HAS  FINISHED  '.18,'  REACTIONS  EXAMINED.'  'HUMBER  FOUND 

0037  6  »  ',16) 

00  38  C 

0039  C  GET  LOGICAL  UNIT  HUMBER  (LU)  OF  DIALOGUE  DEVICE. 

0040  CALL  RMPAR!  IP) 

0041  C 

0042  C  OBTAIN  LU  OF  INPUT  DEVICE  (CURRENTLY  THE  PROGRAM  IS  SET  UP 

0043  C  TO  READ  THE  INPUT  FILE  OFF  A  CASSETTE  THAT  CAN  BE  CHANGED 

0044  C  HITHOUT  TOO  MUCH  DIFFICULTY.) 

0043  URITE!  IP!  1  ).  100  ) 

0046  READ!  I P  *.  1  ),*)  IN 

0047  C  OBTAIN  LU  OF  LIST  DEVICE 

0048  HRITEC  IP!  1  ),  113  ) 

0049  READ!  IP! 1  ),«  )  LP 

0030  C  NOTE  THAT  OUTPUT  DEVICE  CAN  BE  TO  A  HAG  TAPE  (LU  >  8  OR  IS 

0031  C  OH  THIS  SYSTEM)  IF  SO,  A  CONTROL  MORD  INTERNAL  TO  THE 

0032  C  PROGRAM,  IS  NEEDED 

0033  ICHTL  *  1008  ♦  LP 

0034  C  URITE  OUT  TITLE 

0033  URITE(LP, 106  ) 
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0036 

0037 

0038 

0039 

0060 

0061 

0062 

0063 

0064 

0063 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0032 

0083 

0084 

ooas 

0086 

0087 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0093 

0096 

0097 

0098 

0099 

0100 

0101 

0102 

0103 

0144 

0103 

0106 

0107 

0108 

0149 

0110 


0002  RXNGN  9:38  AH  THU..  16  JULY.  1981 


C  READ  IN  0A7A  FROM  CASSETTE 

C  SEE  SAMPLE  DATA  FILE  FDR  FORMAT  . 

RE  AO (  IN.*)  H 
DO  13  I >1  .  H 

READ(  IN  < 1 04  )  ( H  AH  E (  12.  I).  IZ»1,4),C(  I  ) .  H  <  I  )  <  0 (  I  ) ,  HZ (  I  ). EHTL<  1  ) 

♦  EHTRC  I  ) 

HRITE(LP. !45)(MAME<  12.  I >. 12*1. 4). Ct  ]).H(  1  ).0(  I  )  , HZ ( I  ) . E  M  T  L  <  I  ) 

♦  EHTRC  1  ) 

EHTRC  I  >  =  EHTRC  I  171040  .0 

IS  CONTINUE 

C  READ  IH  HUHBER  OF  THE  DIFFERENT  TYPES  OF  BONDS  FROM  DATA  FILE 

C 

HRITECLP. 112) 

DO  204  I  =  1.  H 

REAOC  IN. Ill)  NC.NCHC  HC )<  NCCCNC  ). HCOCNCl.HOQC  HCl.NQHCHC  ) 

240  CONTINUE 

HRITECLP.  141)  C  I.NCHC  I  ).NCCC  I  ).NCOC  I  )  ,  N  0  0  C  I  )  .  NO  H  C  I  ).  1  =  1. N) 

C 

C 

HRITECLP,  108  ) 

C 

C  LOOP  THROUGH  ALL  POSSIBLE  COMBINATIONS  OF  THE  INPUT  SPECIES. 

C 

00  50  1=1 , N 
00  40  J  =  1  ,  N 
DO  30  K  =  1  .  H 
DO  20  L=  l  .  N 

C  CHECK  FOR  HASS  BALANCE 

T1=C  CCK  )  ♦  C  C  L  )  )  -  C  C  C  J  )+C< I  )  ) 

IFCABSC  T1  )  GT  0  .01  )  CO  TO  20 
T2  =  CHCK  )+HCL))”CHCJ  )  ♦ H  < I  )) 

IFC ABSC T2  )  GT  0  01  >  CO  TO  20 
T3=C0CK)+0CL))-CQCJ)+0CI)> 

IFC  ABSC  T3  )  GT  .  0  01  )  GO  TO  20 
C  CHECK  IF  REACTION  HAS  BEEN  DONE  BEFORE 

IFCK.GT.L.OR.  I  .GT . J  )  GO  TO  20 

C  CHECK  FOR  SIMPLE  PERNUTAT  IONS  OF  REACTAHT7PR0DUCT  PAIRS 

IFC I  .  E9  K  .  AND  .  J  ED . L  )  GO  TO  20 

IFC I . EO  L  AND.  J  EO  L  )  GO  TO  20 

IFC l . EO  K . AND . J  EO  .  K  )  GO  TO  20 
IFC I  .  EQ  .L  .  AND  .  J  ED  .  K  )  GO  TO  20 
II«I  HI 

C  DON'T  ALLOU  REACTIONS  OF  STABLE  SPECIES. 

C  THIS  IS  AN  EXAMPLE  OF  A  REASONABLE  CRITERION  THAT  THE  USER 

C  NIGHT  IMPOSE  BASED  ON  PRIOR  KNOHLEDGE  OF  THE  SYSTEM. 

C  NOTE:  HERE  THE  DATA  FILE  HAS  BEEN  SET  UP  IN  A  SPECIAL  HAY 

C  TO  GROUP  ALL  STABLE  SPECIES  SEPARATELY  FROM  ALL  RADICALS 

C  IE.,  THE  RADICALS  HAVE  SPECIE  NUMBERS  >  27. 

IFCHZCI)  GT.  27  AND  HZCJ)  CT  27)  GOTO  20 
C 

IA  *  NZC  I  ) 

JA  •  HZCJ) 

KA  *  NZC  K  ) 

LA  •  HZCL) 
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om 
0112 
0113 
01  14 
0  115 
0  116 
oil? 
0  118 
ol  13 
0  12  0 
0  12! 
0122 
0  123 
0  124 
0  125 
0  126 
0  127 
0  128 
0  129 
0  1  30 
0131 
0132 
0  133 
0134 
0  1  33 
0  136 
0137 
0138 
0  139 
0  140 
0141 
0142 
0  143 
0144 
0145 
0  146 
0147 
0  148 
0149 
0130 
0151 
0132 
0133 
0134 
0135 
0136 
0137 
0158 
0139 
01(0 
01(1 
01(2 
01(3 
0  1(4 
0  1(5 
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C  ALLOW  H  2<  5  3  )  <  C02<40)  AND  H2DM9)  AS  PRODUCTS  ONLY 

!F<  IA  .EQ  55  OR.  IA  EQ.  40  OR.  1A  .  EQ  .  49)  GOTO  20 
IF (  J A  .EQ.  55  OR.  JA  EQ  40  OR.  JA  .  EQ  49)  GOTO  20 
C 

C  FIND  NZtSPECIES  M)  --  SPECIES  KUHBER  OF  THE  SPECIAL  SPECIES  N 

C 

00  201  IFH  •  1,  H 

IF(CCIFH)  HE.  0  .OR.  O(IFH)  HE  0  OR.  H(IFH)  .HE.  0)  GOTO  201 
HM  =  H  2  C IFH) 

GOTO  2o 2 

20 1  COHTIHUE 

202  COHTIHUE 

C  ALLOW  ISOHERIZATIOHS  OHLY  WITH  H 

IFMIH  HE  HH  AMD . JA  NE  HH  I  AND  (  IA  EQ  .KA  OR.  IA.EQ.LA  OR. 
t,  JH.EQ.KA  .OR.  JA.EQ.LA>>  GOTO  20 
C  OIILY  ALLOW  THF<  30  ),  1  -BUTEHEE  36  )  AMD  T-BUTEHEl  35  )  AS  PRODUCTS 

C  OF  tSOHERIZATIONS 

I  r  (  <  K  A  Efl  39  OR  K  H  E  0  3  6  OR.  KA  EQ  33  OR.  LA.E0  38  OR. 
t  LA  E  0  36  OR.  LA.EQ.35)  A  HD  .  (IA.HE.HH  AMD  JA  HE  HH))  G0TO2O 
C 


C 

C 

c 


c 

c 

c 

c 

c 

C 

C 


c 


c 

C 


oei.h  =  i:ehtl<k)*entl<l>)-i'entl<i>  +  ehtl(J)) 

DELS  »  <EHTR(K)*EHTR(L))  -  (EHTR(I)+EHTR(J)> 
DELG  =  DELH  -  TEHP  ♦  DELS 

IFtDELG  .  G  T  20.0  .OR  DELH  .  G  T  .  35.0)  COTO  20 


HAKE  CHECKS 


BY  THE  NUMBER  OF  BONDS  FORHED  AMD  BROKEN 


HCHD  =>  HCH(KA)  ♦  HCH(LA)  -  HCH(IA)  -  HCHCJA) 

NCCO  =  HCC(KA)  >  HCC(LA)  -  HCC(IA)  -  HCC(JA) 

HCOD  =>  II C  0  ( K  A  )  +  HCOCLA)  -  HCQCIA)  -  HCO(JA) 

HOOD  *  H  0  0  <  K  A  )  *  HOOYLA)  -  HOO(!A)  -  HOOfJA) 

HOHO  =  KOfKKA)  *  HOH(LA)  -  NOHCIA)  -  HOIKJA) 

II ET  =  IA8S(  HCCD  )  +  I A  8  S  f  HCHD  ‘•►lABSCMOOD  )  *  I  A  B  S  (  HOHO)*IABS(  HCOD  ) 


SET  SOHE  CRITERIA  FOR  HAXIHUH  HUMBER  OF  BOND  CHANGES  ALLOWED 
IH  A  SPECIFIC  REACTION. 


IFOIET  GT.  3)  COTO  20 


ALLOWS  3  BOND  CHAHGES  IF  CO2<40>  IS  FORKED  OR  IF  2  RADICALS 
FORM  2  STABLE  SPECIES. 

IF< < HET . EQ . 3  ) . AND . < KA . E0 . 40  .  OR  .  L A  .  EQ  .  40  )  )  COTO  90 

IF (<HET. EQ3). A  HD. <IACT. 27. ORJA. GT.27.aR.KA.LT. 28. OR.LA.LT. 28)) 
6  GOTO  20 
90  CONTIHUE 


IF ( IABSCHCOD  ) 

.  CT  . 

2  1 

COTO 

20 

IF< I ABSCHCCD  ) 

GT  . 

1  ) 

COTO 

20 

IF( I  A  8  S  (  M  C  HD  ) 

.  GT  . 

2  ) 

GOTO 

20 

IF<  IABSOIOOD  ) 

.  CT  . 

2  ) 

GOTO 

20 

IF(  I  ABS(  HOltD  ) 

GT  . 

t  ) 

GOTO 

20 

ADD  SPECIFIC  CHECKS 

ALLOW  COM4'  TO  GO  TO  OHLY  C02(  40  )  OR  H<4) 

IF(<IA.E8  44.0R.JA  EQ  44)  AHD.MA  NE.40.AHD  KA.HE.4AHD  LA.HE.4 
&  AMD  LA  HE . 40  )  >  GOTO  20 
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PACE 


0  1  6  6 
0  157 
0  16  0 
0  169 
0  1  70 
0171 
0172 
0173 
0174 
0175 
0176 
0177 
0178 
0179 

oiao 
0181 
0  182 
0183 
0  184 

0  1  as 
0186 
0  187 
0138 
0139 
0  190 
0191 
0192 
0193 
0194 
0193 
0196 
0197 
0198 
0199 
0200 
0201 
0202 
0203 
0204 
0205 
0206 
0207 
0208 
0209 
0210 
0211 
0212 
0213 
0214 
0213 
0216 
0217 
0218 
0219 
0220 


0004  RXNGN  9:38  AH  THU..  16  JULY.  1981 


C  All.  CM  C II  4(30)  TO  GO  TO  ONLY  CH3<26) 

IF<(  1  ft . E8  .  30  OR  JA  EG  50)  AND  ( KA  HE  26  AND. LA  HE  26))  GOTO  20 
C  ALLOW  C2H6C30)  TO  GO  TO  ONLY  C2H5C7) 

IF ( (  1  A  ES  30  .OR  JA  EO  301. AND  ( KA  HE. 7  .AND. LA. HE. 7  ))  GOTO  20 
C  ALLOW  C3 II  0(54)  TO  GO  TO  ONLY  C  H  3  C  H  2  C  H  2  <  1  1  )  OR  CH3CHCH3(12) 

IF((  IA.EQ.54  OR  JA  EG  54)  AND  (KA.HE.ll  A  HD  K A  HE  12  AND. 

&  LA  HE .  1  1  AND  .  LA  .  HE  .  12  )  )  GOTO  20 

C  ALLOW  C4 II  10(37)  TO  GO  TO  ONLY  C II  3C  H  2C  HC  H  3  (  1  5  )  OR  C  H  3  C  H  2  CH  2C  H  2  (  1  4  ) 

IF ( (  I  A  £Q . 37  OR  .  JA , EO  37  )  AND  ( KA  HE  14.AHD.KA  NE.15  AND 
6  LA. HE.  14  AND  .  LA . HE  .  15  )  )  GOTO  20 
C  ALLOW  CH30C5)  TO  GO  TO  SOMETHING  WITH  ZERO  OR  ONE  CARBON 

IF(<  I  A  EO . 5 . OR  .  JA  EQ  5  )  AND  .  (  C(  K )  HE  0  AND . C ( K  )  HE.l  AND. 

&  C ( L )  HE . 0  AND  C(L).NE.I))  GOTO  20 
C  ALLOW  11(1)  TO  GO  TO  ONLY  N(4).  H  2  <  5  5).  H20(  49  >.  OR  H202(  4  1  ) 

!F((IA.EQ.1.QRJA.EQ  1  )  A  NO  <  K  A  HE. 4. AND. KA.  HE. 35. AND. KA. HE. 49 
G  AND  K A  HE  41  AHO  LA. HE. 4  AND  LA  HE  55  AND  LA  HE  49 

&  AND  LA  HE  41  )  )  G  0  T  Q  24 

C  ALLOW  0  H( 2  )  TO  CO  TO  ONLY  H<4>,  H20(49  ).  HQ2(3>,  AND  H202(4  1  ) 

I F ( <  I A  EQ.2  OR  JA  EO  2)  AND  (KA  HE  . 4  AND  KA  HE  49  AND.  K  A  HE  3  AND. 
G  KA  ME  41  AND  LA  HE  4  HUD  LA  HE  49  AMD  LA  HE  3  AMD  LA  HE  41)) 

G  GOTO  20 

C  DOH’T  ALLOU  CERTAIN  CROUPS  OF  ISOMERS  TO  INTERCONNECT  UNLESS  H(4) 

C  IS  INVOLVED  GROUP  1  (  1  1  ,  52.  2  5  ,  22  I  WILL  NOT  CONNECT  WITH 

C  GROUP  II  (12.  53.  13,  23).  GROUP  III  (14,  16,  20,  47)  WILL 

C  NOT  CONNECT  WITH  GROUP  IV  (15,  18,  27-  46). 

1F(( IA.E6  11  OR  IA  EQ.52  OR  IA  EO  25  OR  IA.EQ  22  OR  JA  EO  11  OR. 

G  JA  EQ  52  OR  .  JA  EO  .  25  OR  JA  EQ.22)  AND  (KA  EQ  l 2  OR . KA . EO  .  S3 
6  OR  KA  EQ.  13  OR  .  KA  EQ  .23  OR  LA  EQ  12  OR  LA  EO  53  OR  . LA  EO  13 

G  OR  LA. EO. 23)  ANO.  (IA  NE.HN  OR.  JA  HE  HH  )  )  GOTO  20 

IF((IA£Q.12.0RIA.EQ  530RIA  EQ13.0R.IA.EQ230RJAE8.120R. 

G  JA  EQ  53. OR  JA  EQ.13  OR  JA.EQ  23)  AND.  ( K A  .  E Q  .  1  1  0 R . K A . EO  .  5 2 
G  OR  .KA  .EO  .  25  OR. KA. EQ.220R  LA  EQ.il. OR. LA. EO. 52. OR.  LA. EQ  25 

t  OR. LA. EO. 22)  AND.  (IA.NE.NN  OR.  JA  HE  HH))  GOTO  20 

I F ( (  I A  E8  14. OR  IA.EQ  16  8R.IA  EO  20  OR  IA  EG  47  OR  JA  EQ  14  OR 

G  JA  EQ.16.0R  JA  EG  20  flR.JA  EQ  4 7  )  .  AHO  .  I KA  EQ  1 5  OR  KA  .  EQ  .  1 8 

G  OR  KA  EQ.27  OR  .  KA  EQ  46. OR  LA.EQ150R.LA  EQ  18  OR. 

6  LA  EQ  .  27  OR  .  LA  .  EG  46  )  AND.  (IA.NE.NN  OR.  JA  HE  NH))  GOTO  20 

I F ( <  I  A  EQ  .  15  OR  .  IA . EG  18  OR  IA.EQ  27  OR. IA.EQ. 46  OR  JA  EQ  15. OR. 

G  JA  EQ  18. OR. JA.EQ  27  OR.JA  EQ.46)  AND.  (KA.EQ.14  0R.KA.E8.16 
G  OR  KA  EQ . 20  .  OR  .  KA  EG  47. OR  LA  EQ  14  OR  LA.E8  16  OR  LA, EQ.  20 

6  OR  LAEQ.47)  AHO.  (IA.NE.NN  OR  JA  HE  NM)>  GOTO  20 

C  ALLOW  METHANE  (50)  TO  COHE  FROM  HETHVl  (26)  ONLY 

I F ( ( K  A  EQ . 50  .OR  LA  EO  30)  AND.  (IA  NE  26  AKD.JA  HE  26))  GOTO  20 
J  J  *  J  J  ♦  1 

WRITECLP, 107)(HAME( 12,  I  ),  IZ*  1 ,4),IA,(HAME( 12, J  )  ,  I  Z  *  1 ,4  ) ,  J  A , 

KHAHE(  12,  K  ),  IZ«1, 4),KA,<  NRHE(  IZ.  L).  IZ»1  .  4  >,LA. 

GOELH.OELG,  II, JJ.NCHD.NCCO, NCOD.NOOD.HOHD.NET 
20  CONTINUE 
30  CONTINUE 
40  CONTINUE 
50  CONTINUE 
C 

C  WRITE  OUT  TOTAL  COHBINATIONS  TRIED 

C 

WRITE! LP, 1  10  )  11 


V. 
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022  1  C 
0222 
0223  C 
0224 
0223 


IF  WRITING  TO  TAPE  WRITE  EOF 

IF(  LP  Efl  8  OR  LP  £8  18)  CALL  E  XE  C  (  3,  ICHTL  ) 

INDICATE  END 

WRITE( !P< 1 >-  114  >  1 1  -  J J 

END 
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Appendix  B 
PROGRAM  SELET 


PACE  000  1  FIN  9:46  AM  F  R  I  ,  1?  JULY,  1981 


0001  M  N  i  I. 

0002  PROGRAM  SELET'3,99),  REV  800  710  D  0  RE  H  IHDRITZ 

0043  C  VERSION  000710  FOR  HEWLETT  PACKARD  MINICOMPUTER 

0044  C  PROGRAM  TO  REND  «  SELECTED  DATA  FILE  FOR  SPECIFIED 

0045  II  fi  E  A  C  T  A  N  '  UUR  •'  p  (1  D II  C  T  l  i  N  E  SPECIES  I'Ni.Yi 

004*  C  THIS  ALLOWS  THE  CHEMIST  T-i  1  fi  I.  L  0  H  i  H  £  REHCT14H5  0  F 

0047  C  A  PARTICULAR  SPECIES  IN  A  LAHsE  MECHANISM 

0048  C 

4049  C  LOGICAL  LOC  TELLS  TRUTH  OF  HATCHES 

4010  LOGICAL  LOG 

0011  DIMENSION  I  F < 5  ) , M AMF ( 3  ) ,  I BUFR<  1  3 2  )  -  I  BUF t  1  3 2 > ,  I  D CBi  1  4 4  ) , 

0012  4  ISLCT.8  ).  I  N  A  S ( A  ) 

0013  C  ICR  IS  THE  CARTRIDGE  THAT  THE  DATA  FILE  IS  OH. 

0014  C  HMAT  IS  THE  HU  H  BE  R  OF  SUCCESSFUL  MATCHES 

0415  C  IPL  IS  ASCII  '-f,  USED  TO  LOCATE  POSITION  IN  FILE. 

001*  C  ISLCT  ARE  THE  COLUMNS  IN  THE  DATA  FILE  C 0 R R E S PO HD  I  M G  TO 

4017  C  THE  SPECIES  NI1H8ERS 

0018  DATA  ICR,  I  SC  ,NHAT  ,  I  PL  .  !L/2?,0, 0, 254008- 1  32/ 

0019  DATA  ISLCT/I2, 13,27,28,43,44,58,597 

0020  C 

0021  C  FORMATS 

0022  C 

0023  100  FORMAT'  *  ENTER  DATA  FILE  NAME  (6  CHARACTERS)') 

0024  141  FORMA  T( 3A2  ) 

0025  102  FORMAT'  '  IER  -  '  -  18  ) 

002*  103  FORMAT'  '  WHICH  SPECIES  NUMBER  ARE  YOU  INTERESTED  IH?') 

0027  1  04  FORMAT'  132A1  ) 

0428  10S  FORMAT'  '  '  -132A1  ) 

0029  106  FORMAT'  ' 1REACTI0HS  INVOLVING  SPECIES  NUMBER', 13/ 

0030  4  '  AS  A  PRODUCT'/,'/) 

0031  1  07  FORMAT'  ' 1  REACT  I  UNS  INVOLVING  SPECIES  HUMBER',  13, 

0032  4  'ASA  REACTANT'///  ) 

0433  1  08  FORMAT'  ' 1  REACT  I ONS  INVOLVING  SPECIES  HUMBER',  13,'  AS  A  REACTANT', 

0034  I  'ORA  PRODUCT ' ///  ) 

0035  1  0?  FORMAT'  ’  DO  YOU  NISH  TO  EXAMINE  SPECIES  NUMBER', 13,'  AS  ', 

003*  4  /'A  REACTANT  ONLY?  ■  1'/'  A  PRODUCT  ONLY?  =  2'/ 

0037  4  '  OR  RS  A  PRODUCT  AND  A  REACTANT?  *  3'  ) 

4038  110  FORMAT'//'  SELECTION  PROCESS  FINISHED  A  TOTAL  OF  ',13, 

0039  4  '  SELECTED  '  ) 

0040  111  FORMAT'  10X, 'REACTANTS' ,22X, 'PRODUCTS' ,IBX, ' D  E  L  H ' , 8  X ,  'DEIC'  , 

0041  4  4X,  1  »  T  R  D  '  ,  2  X  ,  '  »  F  M  D  '  ,  3  X  ,  '  H  C  H  0  1  ,  IX,  '  N  C  C  D  '  ,  IX,  '  N  C  0  D  '  ,  IX, 

0042  4  'HOOD' ,  IX,  'NOHO' - 2X,  'NET' // ) 

0043  C 
0044  C 

0 V 4 5  C  GET  LOGICAL  UNIT  NUNBER'LU)  OF  DIALOGUE  DEVICE 

004*  CALL  RNPAR'IP) 

0047  C  SET  LU  FOR  LIST  DEVICE 

0448  LP  *  0* 

0849  C  GET  NAME  OF  DATA  FILE 

0050  URITF' IP' 1  ), 100  ) 

0451  READ'  IP'I  ),  101  )  ( MAMF' K  )  ,  K  >  1,  3> 

0052  C  GET  THE  NUMBER  OF  THE  SPECIES  TO  DO  THE  SEARCH  OH 

0453  WRITE' IP'  1  ),  103  ) 

0054  C  NOTE  *•*  FORMAT  IS  FREE  FORM  INPUT 

4055  READ'  IP'I  ),*  )  I  HATCH 
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0036  C 
0037  C 
0038  1 

0039 
0060  C 
0061 
0062  C 
0063 
0064 
0063 
0066 
0067 
0068  C 
0069  C 
0070 
0071  C 
0072 
0073  C 
0074  C 
0075  7 

0076  C 
0077 
0078  C 
0079 
0080  C 
0081  C 
0082 
0083 
0084  C 
0083  C 
0086  C 
0087 
0088  C 
0089  C 
0090  C 
0091  C 
0092 
0093 
0094 
0093 
0096 

0097  50 

0098  C 
0099 
0140 
0141 
0142 
0143  C 
0144  C 
0143  C 
0146 
0147  C 
0148  C 
0149  C 
0110 


ASK  IF  USER  WISHES  SEARCH  AS  REACTANT  ONLY,  PRODUCT  OH  L  Y  OR 
OR  AS  REACTANTS  AND  PRODUCTS 
BR1TEC 1PC 1), 109  '  IMATf  H 
READC  IPC 1  ) ,  *  )  1 BR 

HAKE  SURE  USER  GAVE  A  LEGITIMATE  RESPONSE 
I F  C  I  8  R  IT  I  OS  IBR  GT  3  1  GOTO  5 

ASSIGN  CORRECT  FORMAT  FOR  OUTPUT  DEPENDING  ON  RESPOHSE 

1 F  *  IBR  EU  1)  ASSIGN  107  TO  ITIT 

I F  C IBR  .  Ed  .  2)  ASSIGN  106  TO  ITIT 

1F< I BR  EO  3)  ASSIGN  108  TO  ITIT 

BR1TECLP- ITIT  )  IMATCN 

URITECLP- 111) 

OPEN  DATA  FILE  ON  DISK  «TCR 

CALL  OPENTIDCB-IER-NAHF.O. ISC  -  ICR- 

CHECK  FOR  ERROR  RETURN 

IF ( I E  R  LT  0)  UR1TEC1PC 11,102)  1ER 

READ  IN  A  LINE  OF  DATA 

CALL  REAOF'  IOCB.IER.  1 8UF R , I L - LEH  ) 

CHECK  FOR  ERROR  RETURN 

1FCIER  LT  0)  URITE<:  IP<  1  )- 102)  IER 

CHECK  FOR  EHD  OF  FILE  -  EOF  ) 

I F ( IER  EQ  -12  )  GOTO  25 

CALL  CODE  READ  DOES  A  FORMAT  TED  READ  OF  THE  FILE. 

CALL  CODE 

READC  IBUFR-104)  <I8UF<I>,  I  «  1-  132) 

LOOK  FOR  FIRST  REACTION  (LOOK  FOR  A  PLUS  SIGN) 

1 F C I AND< IBUFC16  ), 1PL  )  HE  IPL  )  GOTO  7 

CONVERT  SELECTED  PORTIONS  FROM  ASCII  TO  NUMBERS  INTERNAL  TO  THE 
COMPUTER  ISLCT  POINTS  TO  CORRECT  COLUMNS  OF  DATA  FILE 

DO  10  I  -  1,  8 
I M  AS (  I  )  «  IBUFC  I  SLC  T  <  I  )  ) 

IHAS<  I )  ■  IANOC  IMAS(  I  )  .  374008  M400B 
IFdMAS(I)  E0  4  0  6)  IHAS(I)  •  IHASU)  ♦  24B 
I M A S (  I )  •  1 H AS (  I)  -  60S 
CONTINUE 

RECONSTRUCT  NUN8ERS  1A  8  JA  FOR  REACT RHTS.  KA  6  LA  FOR  PRODUCTS 
1A  *  10  *  IHAS(  1  )  ♦  I H AS < 2  ) 

JA  *  10  *  IHAS( 3  )  ♦  IHASC 4  ) 

KA  *  10  *  IHAS<  3  )  ♦  IHASC  6  ) 

LA  •  10  •  1 M  AS ( 7  )  ♦  IHASC  8  ) 

LOOK  FOR  MATCHES 

COTOC II ,12,13),  IBR 

REACTANT  ONLY 


1 

i 

| 


] 


i 


I 


11  LOG  ■  IMATCN  E9  IA  OR  IMATCN  EQ  JA 


INDRITZ,  MADAY,  SHEINSON,  AND  SHAUB 


PAGE 

000  3 

SELET  9:46  AH  PRI  .  17  JULY.  1981 

f 

0)11 

COTO  20 

0112 

C 

0113 

c 

PRODUCTS  ONLY 

0114 

c 

0115 

1  1' 

LUC  »  IMATOH  10  K  A  OR  I M A T C H  EO 

LA 

01  16 

GOTO  20 

0117 

c 

01  16 

c 

REACTANT  670R  PRODUCT 

0119 

c 

0120 

13 

LUC  ’  INATCH  Ed.  1A  OR  1  MAT  CH  Efl 

JA  OR 

INATCH  Efl  KA 

0121 

6  OR  I  HAT  CH  EQ  LA 

« 

0122 

C 

1 

0123 

c 

OUTPUT  TO  LIST  DEVICE 

i 

0124 

c 

01  23 

.!  0 

I F ( L  0  G  )  NRIT£(LP,  103)  (IS OF (I).  I  =  l, 

1  26  ) 

1 

i 

0  126 

I  F  <  LOG  '■  N  K  A  T  *  NMAT  ■»  1 

0127 

C 

| 

01  26 

C 

CHCOK  FOR  EOF 

0129 

C 

0130 

IF'IER  NE  -12)  GOTO  7 

0131 

25 

COHT I  HUE 

0132 

C 

0133 

C 

EOF 

0134 

c 

0133 

c 

HR1TE  HUMBER  OF  MATCHES  TO  LIST  DEVICE 

A  HD  TO 

DIALOGUE  DEVICE. 

0136 

NR  ITEiLP.liO)  UR  A  T 

0137 

M  R  I  T  E  C  I  P  <  1  '  ,  1  1  0  )  H  M  A  T 

0138 

END 

FTM4  COMPILER:  HP 92 0 « 0 -  1  6 09 2  REV  2026  (  800423) 
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Appendix  C 
PROGRAM  SLR  XT 


MCE 


0001 

0042 

0003 

0004 

0005 

0006 

0007 

0008 

0008 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0018 

0020 

0021 

0022 

0023 

0024 

0025 

0026 

0027 

0028 

0028 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 


0001  FTH  8!57  AH  F  R  I  .  ,  17  JULY,  1981 


F  1 N  4  ,  L 

PROGRAM  S  L  R  X  T!  3,99),  REV.  800  8  1  3  DQREH  IHDRIT2 
C  VERSION  800813  FOR  HP  MINICOMPUTER 

C 

C  PROGRAM  TO  READ  A  SELECTED  DATA  FILE  FOR  SPECIFIED 

C  REACTANTS  WAR  PRODUCTS 

C  THIS  PROGRAM  ALLOWS  ONE  TO  SELECT  MULTIPLE  INPUTS  C  <  =  2  0  )  AS 

C  REACTANTS.  THUS  ALLOWING  ONE  TO  SELECT  CLASSES  OF 

C  REACTIONS  (DIFFERS  FROM  GSELET,  WHICH  ALLOWS  ONLY  0  H  E  ) 

C 

C  LOG  IS  LOCICAL  FOR  REACTANT  SPECIES  HATCH 

C  L0G2  IS  LOCICAL  FOR  COMPLEMENTARY  REACTIONS 

LOCICAL  L 0 C , LOG 2 

DIMENSION  IP(5  )  ,  N  A  M  F  <  3  ) ,  IBUFR<  132),  IBUFC  13  2),  IDCB<  14  4), 

6  ISLCT(S),IHAS(8),IIN<20),JIH!2O> 

DATA  ICR,  I  S  C  ,  N  M  A  T  ,  I  P  L  ,  IL/2?,0, 0,254008,  132  /,  IYES/2HYE/ 

DATA  ISLCT/12, 13, 27, 28,43,44, 58, 59/ 

C 

C  FORMATS 

C 

100  FORM A  T (  '  ENTER  DATA  FILE  NAME  (6  CHARACTERS)') 

101  FORMA  T<  3A2  ) 

102  F  ORH A  T (  '  IER  •  '  -  18  ) 

103  FO  RM A  T  <  '  HOW  MANY  ■ I  A ■  REACTANTS  ARE  YOU  INTERESTED  IN?’) 

104  FORMA  T(  132A1  ) 

105  FO  RM A  T<  '  '  ,132A1  ) 

106  FORMAT!  '  ENTER  SPECIES  ->  ') 

107  FORMAT!  ' 1REACTI0NS  INVOLVING  SPECIES  AS  REACTANTS-)  ',  ,  1 X , 2  0  13) 

108  FORMAT!'  00  YOU  WANT  TO  SPECIFY  COMPLEMENTARY  REACTANTS?') 

109  FORMAT! A2  ) 

110  FORMAT!//'  SELECTION  PROCESS  FINISHED.  A  TOTAL  OF  ',13, 

6  '  SELECTED  '  ) 

111  FORMAT!  10K, 'REACTANTS'  ,22N,  'PRODUCTS',  1  8  M ,  '  D  E  L  H ' , 8X  ,  '  D  E  L  G  ’  , 

6  4X,'»TRD',2X,'»FND',3X,'NCHD',1X,'NCCD'.1X,'HC0D',1X, 

6  'HOOD' , IX, 'NOHD' , 2X, 'NET' // ) 

112  FORMAT!'  THIS  PROGRAM  ALLOWS  YOU  TO  SELECT  UP  TO  20  REACTANTS') 

113  FORMAT!'  HOW  MAMY  -JA*  REACTANTS  ARE  YOU  INTERESTED  IN-*') 

114  FORMAT!'  COMPLEMENTARY  REACTANTS  ->  ',2013) 

C 

C 

C  OBTAIN  LOGICAL  UNIT  HUN8ER! LU  )  OF  INPUT  OEVICE  FOR  DIALOCUE 

CALL  RMPAR! IP) 

C  LP  IS  LU  FOR  OUTPUT  LISTING 

LP  •  06 

WRITE!  IP!  1  ),  112  ) 

C  OBTAIN  NAME  OF  DATA  FILE 

WRITE!  IP!  1  ), 100  ) 

READ!  IP!  1  ) ,  1  01  )  ! HAHF! K  >  ,  K  >  1,  3) 

C  OBTAIN  HUMBER  OF  REACTANTS 

WRITE!  IP!  1),  103  ) 

READ!  IP! 1  ),*  )  Nil 

C  OBTAIN  NUMBERS  OF  THOSE  REACTANTS 

C  NAVE  TO  REFER  TO  ORIGINAL  DATA  INPUT  FILE  TO  SEE  SPECIES  NUMBERING 

WRITE! IP! 1  ),  106  ) 

READ!  IP! 1  ) - •  )  i  1 1  H! L  L  ),  LL  ■  1,  Nil) 
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003«  C 
0057 
0039 
0039 
0060 
00(1  C 
0  0(2 
00(3 
00(4  C 
00(5 
00(( 

00(7  9 

00(8  C 
00(9 
0070 
0071 
0072  C 
0073  C 
0074 
0075  C 
0076 
0077  C 
0078  C 
0079  7 

0080  C 
0081 
0082  C 
0083 
0084  C 
0085  C 
0086  C 
0087  C 
0088  C 
0089 
0090 
0091  C 
0092  C 
0093  C 
0094 
0093  C 
0096  C 
0097  C 
0098  C 
0099 
0  100  C 
0101  C 
0102 
0103 
0104 
0103 

0106  10 
0107 
0108 
0109  C 
0110  C 


SEE  IF  USER  BAHTS  TO  SPECIFY  COMPLEMENTARY  REACTANTS 
BR ITE< IP<  1  >,  100  > 

REACT IP<  1  ) •  1  09  )  I  A H  S 
I F ( IANS  E8  IVES  )  L0C2  =  TRUE 
I F (  NOT  L0C2  )  GOTO  9 

OBTAIN  NUH8EA  OF  CO H P L EH E N T AR Y  REACTANTS 
BR I  TEC IP< 1  I, 113  ) 

REACT  1PT l  ),»  )  HJJ 

OBTAIN  THE  NUH6ERS  OF  THOSE  REACTANTS 
BRITET IPT  1  ),  106  ) 

REACT  1PT  1  )  ,*  )  IJINLLL).  LI  •  1,  HJJ) 

CONTINUE 

BRITE  HEADER  UH1CH  INCLUDES  REACTANT  SPECIES  NUMBERS 
BRITE' LP.  107  1  T  I  I  NT L  ), L  =  1  .Nil  ) 

IFTL0G2)  BRITET LP -  1  14  )  LJ1NLL).  L  =  I  ,  HJJ) 

URITETLP, 111) 

OPEN  OATA  FILE  FROM  DISK 

CALL  OPENT  IDCB>  IER/NAHF. 0 • I  SC  «  ICR) 

CHECK  FOR  ERROR  RETURN 

IFTIER  LT  0)  BRITE)  IP'  1  ),  102 >  I E R 

READ  IN  DATA  LINE  AT  A  TINE 
CALL  AEADFT1OC0,IER,IBUFR,!L,LEN) 

CHECK  FOR  ERROR  RETURN 

IFTIER  LT  0)  BRITET  IPT  1  ),  102  )  IER 

CHECK  FOR  END  OF  FILE  TEOF) 

IFTIER  E  8  -12  )  GOTO  25 

READ  FILE  BITH  FORMAT  NOTE  THAT  A  READF  CALL  GETS  RECORD 
FROM  FILE  AS  IS  CALL  CODE  READ  DOES  A  FORMATTED  READ  OF  THE 
FILE.  HP  IS  SUPPOSED  TO  CONE  OUT  BITH  A  SOFTBARE  REVISION 
THAT  ALLOBS  FOR  DIRECT  REAOING  OF  A  NAMED  FILE 
CALL  CODE 

READ! IBUFR , 104  )  T  I B  U  F  T  I  )  <  I  >  1,  132) 

LOOK  FOR  FIRST  REACTION  T  LOOK  FOR  A  PLUS  SIGH) 

IFT IANDTIBUFTI6  ),  IPL  )  HE  IPL)  GOTO  7 

CONVERT  SELECTED  PORTIONS  FROM  KNOHK  FORMAT  INTO  NUMBERS 
INTERNAL  TO  THE  COMPUTER 

00  10  I  -  1,  8 

FOLLBIHG  LINE  TOGETHER  BITH  DATA  STATEMENT  FOR  ISLCT  GETS 
CHARACTERS  FROH  SPECIFIC  COLUMNS  OF  DATA  FILE 
1 H  AS (  1 )  -  I  B  UF (  ISLCT  <  I  )  ) 

IMAS(I)  •  IAHD<  IHASL  I  )  ,374008  17400B 
IFLIHASLI)  EO  408  )  IMAS(l)  -  IHAS(I)  ♦  20B 
1 M  A S (  I  )  ■  I M  AS  <  1  )  -  (OB 
CONTINUE 

IA  •  10  •  IHASL  1  >  ♦  IHASL  2  ) 

JA  »  10  *  IHASL 3  )  ♦  IHASL  4  ) 

KA  ANO  LA  ARE  THE  NUMBERS  OF  THE  PRODUCTS  FOR  EACH  REACTION 
THEY  ARE  NOT  USED  IN  THIS  PROGRAM  THE  INFORMATION  IS  INCLUDED 


24 


NRL  REPORT  8498 


PAGE 


01 11 
0112 
0113 
0114 
0113 
01  It 
0117 
0118 
0119 
0120 
0121 
0122 
0123 
0124 
0123 
0126 
0127 
0128 
0129 
0U0 
0U1 
0  132 
0133 
0134 
0133 
01  3« 
0137 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0143 
0  14  6 
0147 
0148 
0149 
0130 
0131 
0132 
0133 
0134 
0133 
0136 
0137 


FTH4 


«• 


0003  »LR*T  9 1 3  7  AH  PM.,  17  JULY.  1981 


C  HERE  SO  THAT  ONE  CRN  EASILY  SEE  HOU  TO  IHPLEHENT  A  SEARCH  FOR 

C  SPECIFIC  PROOUCTS . 

C  KA  =  10  •  IHAS<3)  ♦  IHAS(6> 

C  LA  “  10  *  I M  A  S  (  7  )  ♦  IHAS<  8  ) 

C 

C 

C  IF  HE  ARE  LOOKING  FOR  REACTIONS  OF  SPECIES  WITH  SPECIFIED 

C  COHPLENEN  TARY  REACTANTS-  JUHP  AHEAD  TO  LINE  21. 

I F  L  0  G  2  1  GOTO  21 
00  IS  II  •  1  -  Nil 
LOG  -  FALSE  . 

C  SEE  IF  SPECIFIED  REACTANTS  OCCUR  IN  POSITION  1A  OR  JA. 

IF  < I  I N<  1 1  )  HE.  IA  AND.  IIH(II)  HE  JA  )  GOTO  IS 
LOG  -  TRUE 
GOTO  20 
IS  CONTINUE 
COTO  22 
C 

C  LOOK  FOR  NATCNES  OF  BOTH  SPECIES  AND  CON  PL  EH E NT AR Y  REACTAHTS 

21  CONTINUE 

DO  30  II  a  l  ,  Nil 
DO  30  JJ  3  i,  NJJ 
LOG  -  FALSE. 

IF(<  I  I H ( I  I  )  EQ  IA  AND.  JIN(JJ)  E8  .  JA)  OR  (IIN(II)  EO.  JA 
4  AND.  JIN(JJ)  EO  IA))  LOC  3  TRUE 
IF(LOG)  GOTO  20 
30  CONTINUE 
C 

C  OUTPUT  REACTION  HATCHED 

C 

20  IF(LOG)  URITE< LP- 105  )  (IBUFCI),  I  >  1,  126) 

C  COUNT  NUMBER  OF  MATCHES 

IF(LOC)  HNAT  -  HN AT  *  1 
C 

C  CHECK  FOR  EOF.  IF  NOT  EOF.  CONTINUE  READING 

C 

22  CONTINUE 

IFEIER  HE  .  -12  )  COTO  7 
25  COHTIHUE 

C 

C  EOF 

C 

C  IF  EOF  HRITE  NUMBER  OF  HATCHES  TO  LISTING  DEVICE  AND  TO  DIALOGUE 

C  DEVICE 

HRITE(LP.llO)  NNAT 
NRITEC IP<  1  ), 110  )  NHAT 
END 


COMPILER!  HP92060- 1 6492  REV.  2026  (800423) 


NO  HANNINGS  ««  NO  ERRORS  «•  PROGRAM  ■  01337  COHNON  «  00000 


k 
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0001 
0002 
0003 
0004 
0003 
0006 
0007 
0008 
00  08 
00  10 
0011 
00  12 
0013 
0014 
0015 
0016 
0017 
0018 
0019 
0020 
0021 
0022 
0023 
0024 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
0040 
0041 
0042 
0043 
0044 
0045 
0046 
0047 
0048 
0049 
0050 
0091 
0092 
0093 
0094 
0095 


0081  FTN  9112  OH  WED  i  22  JULY,  1981 


FTN4,  L 

PROGRAM  RLIST!3,99>,  REV  810420  DOREH  IHDRIT2 
C  VERSION  810420  FOR  HEWLETT  PACKARD  MINICOMPUTER 

C  TO  LIST  REaCTIOH  FILE.  DOUBLE  SPACES. 

C  IN  SHORTENED  FORMAT  (ONLY  FIRST  60  COLUMNS  PLUS  THE 

C  NUMBER  LISTED  UNDER  *  AF  H  D  *  )  THIS  GIVES  A  LISTING 

C  THAT  HAS  ROOM  TO  HAKE  COMMENTS  IN  THE  RIGHT  HAND  ’MARGIN*. 

C  BUT  STILL  KEEPS  REQUISITE  IDENTIFIERS 

C  ECHOS  DATA  AT  HEAD  OF  FILE  UNTIL  FIHDS  FIRST  REACTION 

C  PROGRAM  COUNTS  NUMBER  OF  REACTIONS  FOUND  IN  THE  FILE. 

C 

DIMENSION  IP<5).NAMF(3  ) .  IBUFR!  132).  I B  UF  <  132).  1 D  CB  <  144).  IM  AS  <  5  ) 
DATA  ICR.  ISC.NHAT.IPL.  IL.  IBLBL/27.0.0.25404B.132.2H  / 

C  ICR  IS  TIE  CARTRIDGE  THAT  THE  DATA  FILE  IS  OH 

C  IPL  IS  ASCII  FOR  *♦*  IH  THE  LEFT  BYTE  . 

C 

C  FORMATS 

C 

100  FORMAT!'  ENTER  OATA  FILE  MARE  (6  CHARACTERS)') 

101  FORMAT!  3  A  2  i 

102  FO  RM A  T  <  '  I ER  »  '.181 

1 03  FORMA  T  <  132A1  ) 

t  0  4  FORMAT!  1H0.60A1  ,  1  X  ,  14) 

105  FORMAT!'  FIHISHEO  A  TOTAL  OF  ',13,'  WRITTEN.') 

C 

C  OBTAIN  LOGICAL  UNIT  NUMBER! LU  )  OF  DIALOGUE  DEVICE 

CALL  RHPARCIP) 

C  SET  LU  FOR  LIST  DEVICE 

LP  ■  9 

C  OBTAIN  NAME  OF  DATA  FILE 

WRITE!  IP!  1  ).  100  ) 

READ!  IP! t  ),  101)  ! MAMF! K  ) ,  K  ■  1,  3) 

C 

C  OPEN  OATA  FILE  FOR  REROIHG. 

CALL  OPEN! IDC8, IER.NAHF.O.ISC,  ICR) 

IF! I E  R  LT.  0)  WRITE! IP! I).  102  )  IER 
7  00  10  I  -  1,  132 
IBUFR!  I  )  -  IRLBL 
10  CONTINUE 

CALL  READFTIDCB.IER,  IBUFR, IL.LEN) 

IF! IER  LT.  0)  WRITE! IP! 1  ), 102  )  IER 
IF! IER  EQ  -12  )  GOTO  25 
C 

C  CALL  CODE  READ  00E8  A  FORMATTED  READ  OF  A  SPECIFIED  BUFFER 

CRLL  CODE 

READ!  IBUFR  .  103  )  !I8UF!I),  I  ■  1.  132) 

C 

C  LOOK  FOR  FIRST  REACTION  !*  *  *  SIGH) 

C 

IF! I  AND! IBWF! 16  ),  IPL  )  EB.  IPL)  GOTO  8 
WRITE! LP, 103  )  !  1BUF! L  ) , L- 1 , 1 32  ) 

GOTO  7 
C 

C  GET  OFM0»  THI8  WILL  IE  PART  OF  THE  OUTPUT 
C 
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PACE 

0002 

RLIST  9  >  1  2  RH  WED  .  -  22  JULV,  1981 

0038 

8 

80  20  I  >  1.  S 

0037 

IHRS(I)  -  IRHD<  IBUF<  85*1  )  <374008  1/4008 

0038 

IF(IHRSd)  EO  408  )  I  HAS  <  t  )  •  I  HAS  <  I  )  4  2oB 

0039 

IHRS (  I  )  -  1 M  A  S  <  I)  -  (OB 

00(0 

20 

COHTIHUE 

00(1 

REND  »  10000*IHAS<1>  +  1  000*1  MAS<2)M0O»IHAS(3>*10*IHAS<4)+lHAS(5> 

00(2 

C 

00(3 

8R1TECLP-  104  )  (lBUF(I)d  -  1.  (0),  HERO 

00(4 

C 

KEEP  TRACK  OF  THE  NUMBER  OF  HATCHES 

00(3 

RHAT  •  MHAT  ♦  1 

00(( 

C 

00(7 

C 

CHECK  FOR  EOF 

00(8 

IF (  I  E R  HE  -12 i  GOTO  7 

00(9 

25 

COHTIHUE 

0070 

C 

0071 

WRITEdP,  105  )  RHAT 

0072 

HRITEC IP< 1  ), 105 )  RHAT 

0073 

EHD 

FTM4  COMPILER:  HP  92 0 « 0  -  1 6 0 9 2  REV  202i  (  800423  ) 

**  HO  WARB1HGS  •*  HO  ERRORS  *•  PR  OCR  RH  »  0081  3  COHHON  *  00000 
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END 


DATE 

FILMED 


:i 


DTIC 


AlX. 


